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The alternating and non-alternating harmonic sums and other algebraic objects of the same equivalence class are 
connected by algebraic relations which are induced by the product of these quantities and which depend on their 
index calss rather than on their value. We show how to find a basis of the associated algebra. The length of the 
basis I is found to be < 1 /d, where d is the depth of the sums considered and is given by the 2nd Witt formula. 
It can be also determined counting the Lyndon words of the respective index set. The relations derived can be 
used to simplify results of higher order calculations in QED and QCD. 



1. INTRODUCTION 

Multiple nested alternating and non-alternating 
harmonic sums S ail ,,, ian (N) [1-3] emerge in per- 
turbative higher order calculations within QED 
and QCD for massless fermions, 
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Here, au are positive or negative integers and 
TV is a positive even or odd integer depend- 
ing on the observable under consideration. One 
calls n the depth and Y^k=i I * I the weight of 
a harmonic sum. Harmonic sums are associ- 
ated to Mellin transforms of real functions or 
SCHWARTZ-distributions f(x) e S'[0, 1] [4] 
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which emerge in field theoretic calculations. Fi- 
nite harmonic sums are related to harmonic 
polylogarithms if&i, ...,&„ (x) [5]. Their 1/(1 ± 
x)-weighted Mellin transform yields harmonic 
sums. The inverse Mellin transform relates 
the harmonic sums to functions of Nielsen in- 
tegrals [6] of the variable x at least for all sums 
of weight w < 4 as shown in [3], and associated 
generalizations for higher weight. Nielsen inte- 
grals are a generalization of the usual polylog- 



arithms [7]. In the limit N — > oo the conver- 
gent multiple harmonic sums, i.e. those where 
ai 7^ 1, yield (multiple) Zeta-values C 0lv .. l0n , 
which are also called Euler-Zagier sums [8] . A 
generalization of both harmonic polylogarithms 
and the Euler-Zagier sums are the nested Z- 
sums [9], which form a Hopf algebra [10,11] and 
are related to Goncharov's multiple polyloga- 
rithms [12]. For a recent review see [13]. 

Higher order calculations in massless filed the- 
ories are either performed in Mellin-JV space 
referring to harmonic sums or in the space of the 
momentum fractions x representing the results in 
terms of NiELSEN-type integrals. The principal 
complexity is determined by the amount of pos- 
sible terms contributing. In the case of the 2- 
loop coefficient functions in momentum fraction 
space [14] 77 different functions occurred, cf. [3]. 
This number compares to the amount of all pos- 
sible different nested harmonic sums up to weight 
w = 4, 80 = 3 W — 1. For the 3-loop anomalous 
dimensions [15] one expects the contribution of 
a wide class of the w = 5 harmonic sums and 
for the 3-loop coefficient functions of the w = 6 
harmonic sums, which means 242 or 728 sums, re- 
spectively. These sums are not independent but 
connected by different kind of relations. In the 
present paper we summarize a first class of re- 
lations recently being discussed in Ref. [16], the 
so-called algebraic relations. 1 It turns out that 
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these relations emerge from the index-structure 
and the multiplication relation of the objects con- 
sidered and are widely independent of other prop- 
erties of the harmonic sums. In this way an equiv- 
alence class of even more objects is defined hav- 
ing the same properties or can be found as special 
cases thereof. One example is the set of the har- 
monic polylogarithms [5]. 

To obtain manageable expressions it is of im- 
portance to apply all these relations through 
which the number of basic functions to be referred 
to is considerably reduced. Experience shows 
that the Mellin space representation yields sim- 
pler expressions in general [18], which would not 
be seen easily working in x space. To obtain as 
simple as possible expressions it is of special im- 
portance because of the fact that data-analyses 
require compact results in x-space, which can be 
obtained using analytic continuations for the ba- 
sic sums [19] and performing a single numeric 
Mellin inversion [20] for the whole problem. 
Since the evaluation of precise analytic continu- 
ations needs special effort any possible reduction 
carried out before is of help. 

2. ALGEBRAIC RELATIONS 

The product of two finite harmonic sums (1) 
yields 
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We introduce the shuffle product LU of a single 
and a general finite harmonic sum 

S ai (N)LlJS bli ..., bm (N) = S aiM ,... ibm (N) 
+s bl (N) + ... + S buba ,... tbm , ai (N) 

(4) 

which is a linear combination of the sums of depth 
m + 1 . The shuffle product of two harmonic sums 



of depth n and m, S aii ..., an (N) and Sf )li ... i b m (N), 
is the sum of all harmonic sums of depth m + n 
in the index set of which a% occurs left of cij for 
i < j and likewise for bk and bi for k < I. As 
an example the shuffle product of two threefold 
harmonic sums is given by 
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Finally one establishes a system of linear equa- 
tions in which the linear elements of the shuffle 
products form the variables and a polynomial out 
of harmonic sums of lower depth forms the in- 
homogeneity. We furthermore consider all index 
permutations. This system contains all algebraic 
relations. In Ref. [16] all solutions for harmonic 
sums up to depth d = 6 were given. This complies 
to the level of sophistication needed to reduce the 
corresponding relations which emerge for mass- 
less 3-loop coefficient functions. 

3. NUMBER OF ALGEBRAICALLY IN- 
DEPENDENT HARMONIC SUMS 

Let us consider the index set of a harmonic sum 
of depth d. One may consider this set as a word 
wora non-commutative product of letters of an 
ordered alphabet 21 = {a, b, c, d, . . .}. Any word 
can be decomposed into three parts 

w = pxs , (6) 

a prefix p, a suffix s, and the remainder part 
x. Among all words w the Lyndon words, cf. 
e.g. [21], are those being smaller than any of its 
suffixes. 

According to a Theorem by Radford [22] the 
shuffle algebra discussed above is freely generated 
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by the Lyndon words, i.e. the length of its basis 
is given by the number of Lyndon words. We 
would like to count the number of Lyndon words 
for index sets, where the same letters can emerge 
repeatedly. The corresponding relation is due to 
Witt [23] and will be called 2nd Witt formula 
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with n = J2k=i nk - Here /i(d) denotes the 
Mobius function. One may derive l n (ni, . . . , l q ) 
using the generating functional 

oo / i \ 

Note that (7) is related to the Gauss-Witt re- 
lation mentioned by Hoffman [24] for the num- 
ber of basic multiple Zeta-values of weight w for 
Vrii > if all cases for fixed weight are summed 
over. An even more strict relation in the inclu- 
sive case has been conjectured by Zagier [8] and 
Broadhurst and Kreimer [25] in the case of 
multiple Zeta- values and verified up to w = 12. 

Let us come back to Eq. (7). We can draw some 
immediate conclusions out of this relation. If the 
numbers rii have no common divisor larger than 1, 
the number of the basis elements compared to the 
number of all objects equals 1/d, where d denotes 
the depth of the index set. In case of common 
divisors larger than 1 we checked that the basis 
is always shorter for all depths up to d = 10, 
see [16]. 



Weight 


# Sums 


# Basic Sums 


Fraction 


1 


2 





0.0 


2 


8 


1 


0.1250 


3 


26 


7 


0.2692 


4 


80 


23 


0.2875 


5 


242 


69 


0.2851 


6 


728 


183 


0.2513 



4. CONCLUSIONS 



The product of finite alternating or non- 
alternating harmonic sums is given by the shuf- 
fle product of harmonic sums and polynomials of 



harmonic sums of lower depth. These representa- 
tions imply algebraic relations between the har- 
monic sums. If one considers all harmonic sums 
associated to an index set {a\, . . . , a^} one may 
express these sums by a number of basic sums. It 
turns out that this number is given by the 2nd 
Witt formula which counts the number of Lyn- 
don words corresponding to the respective index 
set. The set of these Lyndon words generates in 
this sense all harmonic sums of this class freely. 
By solving the corresponding linear equations we 
derived the explicit representation of all harmonic 
sums up to depth d = 6 without specifying the 
indices numerically and gave all expression which 
are structurally needed to express the sums up to 
weight w = 6. The counting relations for the ba- 
sis of the finite harmonic sums were given up to 
depth d = 10. The relations derived hold likewise 
for other mathematical objects obeying the same 
multiplication relation or a simpler one, which is 
being contained, as that for harmonic polyloga- 
rithms. This is due to the fact that the relations 
derived depend on the index set and the multi- 
plication relation but on no further properties of 
the objects considered. 

The ratio of the number of basic sums for a 
given index set and the number of all sums is 
mainly determined by the depth d rather than 
the weight of the respective sums, due to the pre- 
factor 1/d in the Witt formula. Modifications 
occur due to common non-trivial divisors of the 
numbers of individual indices in the set being con- 
sidered. Up to d = 10 we showed that the frac- 
tion of basic sums is always < 1/d compared to all 
sums. The use of these algebraic relations leads 
to a considerable reduction in the set of functions 
needed to express the results of higher order cal- 
culations in massless QED and QCD and related 
subjects. For practical applications such as the 
description of the QCD scaling violation of the 
structure functions in deeply inelastic scattering 
the harmonic sums occurring in the Mellin N 
space calculation have to be translated to x-space 
by the inverse Mellin transform. For this reason 
the respective harmonic sums have to be analyt- 
ically continued in the argument N to complex 
values, which requires a high effort using numer- 
ical procedures. It is therefore recommended to 
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use as many as possible relations between the N 
space objects before to perform the last step only 
for the reduced set. 

Acknowledgment. This paper was supported 
in part by DFG Sonderforschungsbereich Transre- 
gio 9, Computergestiitzte Theoretische Physik. 

REFERENCES 

1. A. Gonzalez- Arroyo, C. Lopez, and F.J. Yn- 
durain, Nucl. Phys. B153 (1979) 161; 

A. Gonzalez-Arroyo and C. Lopez, Nucl. 
Phys. B166 (1980) 429. 

2. J. A.M. Vermaseren, Int. J. Mod. Phys. A14 
(1999) 2037. 

3. J. Bliimlcin and S. Kurth, DESY 97-160, 
hep-ph/9708388; Phys. Rev. D60 (1999) 
014018. 

4. L. Schwartz, Theorie des Distributions (Her- 
mann and Cie, Paris, 1951), Vols. I and II; 
V.S. Vladimirov, Gleichungen der Mathematis- 
chen Physik, (DVW, Berlin, 1972); 

K. Yoshida, Functional Analysis (Springer, 
Berlin, 1978). 

5. E. Remiddi and J. Vermaseren, Int. J. Mod. 
Phys. A15 (2000) 725. 

6. N. Nielsen, Nova Acta Leopold. XC, 121 
(1909); 

S. Kolbig, Siam J. Math. Anal. 17 (1986) 
1232. 

7. L. Lewin, Dilogarithms and Associated Func- 
tions (Mac-Donald, London, 1958); Polyloga- 
rithms and Associated Functions, (North Hol- 
land, New York, 1981); 

A. Dcvoto and D.W. Duke, Riv. Nuovo Cim. 
7 N6 (1984) 1. 

8. L. Euler, Novi Comm. Acad. Sci. Petropol. 
20 (1775) 140; 

D. Zagier, First European Congress of Math- 
ematics, (Birkhauser, Boston, 1994), Vol.11, 
p. 497. 

9. S. Moch, P. Uwcr, and S. Weinzierl, J. Math. 
Phys. 43 (2002) 3363; 

S. Weinzierl, Comput. Phys. Commun. 145 
(2002) 357. 

10. H. Hopf, Ann. of Math. 42 (1941) 22; 

J. Milncr and J. Moore, Ann. of Math. 81 
(1965) 211; 



M.E. Sweedler, Hopf Algebras, (Benjamin, 
New York, 1969). 

11. C. Kassel, Quantum Groups, (Springer, Berlin, 
1995). 

12. A.B. Goncharov, Math. Res. Lett. 5 (1998) 
497, 

http : // www . math . uiuc . edu/K-Theory/0297. 

13. M. Waldschmidt, J. Theor. Nombres Bor- 
deaux 12 (2000) 581; Multiple polyloga- 
rithms, Lectures at the Institute of Math. Sci- 
ences (Chennai, Nov. 2000); updated version 
http://www.matli.jussieu.fr/~iniw/ 
articles/ps/mpl .ps; 

V.V. Zudilin, Uspckhi Mat. Nauk (Russian 
Math. Surveys) 58 (2003) 3. 

14. W.L van Neerven and E.B. Zijlstra, Nucl. 
Phys. B382 (1992) 11. 

15. S.A. Larin, T. van Ritbergen, and J. A.M. 
Vermaseren, Nucl. Phys. B427 (1994) 41; 
S.A. Larin, P. Nogucira, T. van Ritbergen, 
and J. A.M. Vermaseren Nucl. Phys. B492 

(1997) 338; 

A. Retey and J. A.M. Vermaseren, Nucl. Phys. 
B604 (2001) 281. 

16. J. Bliimlcin, Comput. Phys. Commun. 159 
(2004) 19. 

17. J. Bliimlcin, in preparation. 

18. J. Bliimlcin and S. Moch, in preparation. 

19. J. Bliimlcin, Comput. Phys. Commun. 133 
(2000) 76. 

20. J. Bliimlein and A. Vogt, Phys. Rev. D57 

(1998) 1; D58 (1998) 014020. 

21. C. Reutenauer, Free Lie Algebras, (Calcndron 
Press, Oxford, 1993). 

22. D.E. Radford, J. Algebra, 58 (1979) 432. 

23. E. Witt, Journ. Reine und Angew. Mathe- 
matik, 177 (1937) 152; Math. Zeitschr. 64 
(1956) 195. 

24. M.E. Hoffman, J. Algebra, 194 (1997) 477. 

25. D.J. Broadhurst and D. Kreimer, Phys. Lett. 
B393 (1997) 403. 



